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Abstract. Although the inflationary paradigm is the most widely accepted explanation for the current 
cosmological observations, it does not necessarily correspond to what actually happened in 
the early stages of our Universe. To decide on this issue, two paths can be followed: first, all 
the possible predictions it makes must be derived thoroughly and compared with available 
data, and second, all the imaginable alternatives must be ruled out. Leaving the first task to 
all other contributors of this volume, we concentrate here on the second option, focusing on 
the bouncing alternatives and their consequences. 

Quoique le paradigme inflationaire soit maintenant communement accepte comme represen- 
tant la meilleure explication des donnees cosmologiques, il n'est pas pour autant possible 
de dire qu’une telle phase soit averee. Pour s’approcher d’une telle conclusion, on peut 
suivre deux chemins differents : on peut explorer les consequences de T inflation pour la 
pousser dans ses derniers retranchements, ou bien, au contraire, etudier en detail les alter¬ 
natives possibles. La premiere option faisant l’objet de la plupart des contributions de ce 
volume, nous nous concentrons ici sur la seconde, et presentons les modeles dans lesquels 
une phase de contraction est suivie d'un rebond conduisant a notre epoque d’expansion. 
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1. Introduction 

Starting out in a dense state some 13.8 billion years 
ago, our Universe and its evolution since this ini¬ 
tial time are well understood, with an initially al¬ 
most scale-invariant, but not quite, spectrum of pri¬ 
mordial perturbations condensing into the presently 
observed large scale structures by means of grav¬ 
itational collapse. The very high densities of the 
early stages provide initial conditions to explain the 
relative amounts of different nuclei, and the ensu¬ 
ing phases, being controlled by well-known physical 
mechanisms, permit to reconstruct, from the cosmic 
microwave background (CMB) observations, the prop¬ 
erties of the last scattering surface. We have arrived 
at the point [1, 2] where cosmological data can be 
used to probe the earliest conceivable phases. 

The most widely accepted paradigm for describ¬ 
ing the earliest phases of the Universe, when the en¬ 
ergy density was a mere few orders of magnitude be¬ 
low the Planck scale, is inflation [3,4], Easily imple¬ 
mented by means of a scalar field, this almost expo¬ 
nentially expanding era rapidly leads to a flat Fried¬ 
mann Lemaitre (FL) spacetime with a very slightly 
reddish spectrum of initial perturbations from which 
the rest of the history of the Universe ensues. As 
is, such a scenario is compatible with all currently 
available data. 

This contribution reviews some properties of some 
non inflationary bouncing models. The first natural 
question that comes to mind before going any fur¬ 
ther is: why should we bother with possible alterna¬ 
tives to a working scenario? There are in fact many 
reasons, the first of which being that the phase of 
inflation is silent relative to the primordial singular¬ 
ity, as we discuss in Sec. 2. below. The second is 
that there is no way we will ever be able to assert 
that a phase of inflation did actually take place ex¬ 
cept through its presently observable consequences. 
But then the question arises as to whether other com¬ 
peting theories could induce similar consequences. 
Thuus, examining all plausible scenarios in detail 
seems to be the only way to assert whether infla¬ 
tion is the unique possibility leading to our observ¬ 
able Universe. In the end, ruling out alternatives, or 
not, increases or decreases our level of confidence 
in inflation until it becomes, if ever, recognized as 
valid beyond any reasonable doubt. As we shall see 


in Sec. 3., there are bouncing alternative explana¬ 
tions to the standard model puzzles of homogeneity, 
flatness, isotropy, horizon and the overproduction of 
relics, as well as many models, some of which are 
listed in Sec. 4., in which those bounces can be im¬ 
plemented. 

Getting a background-compatible model is how¬ 
ever not the end of the story: the recently released 
Planck data [5, 6] confirm what was suggested by 
previous experiments, namely that the spectrum of 
primordial perturbations was almost scale invariant: 
slightly red, with a spectral index n s = 0.9639 ± 
0.0047, excluding exact scale invariance at the 5er 
level. The level of non-gaussianity is compatible 
with zero, and the contribution of tensor modes re¬ 
mains below the ~ 10% limit relative to the scalar 
amplitude. All these facts are compatible with the 
perturbations having been produced by quantum vac¬ 
uum fluctuations of a single scalar degree of free¬ 
dom, a natural consequence of slow-roll single-field 
inflation. Can a non-inflationary bouncing repro¬ 
duce such results? As of now, there is no definite 
answer to this question. For this reason, and for lack 
of space in the present article, we shall not discuss 
these points below, and instead refer the reader to a 
recent review [7] in which all the relevant constraints 
for the models exhibited below are derived. 


2. The singularity 

The fact that cosmology, or at least its classical 
implementation in terms of general relativity (GR), 
always leads to the existence of singularities stems 
from the well-known singularity theorems [8], A 
general argument was proposed in Ref. [9]: in an 
FF spacetime with metric 


ds 2 = —dt 2 + a 2 (t) 7^ (a;) dx'daP 

= a 2 (?y) [—d?7 2 + 7^ (a;) dx z dx j ] , (1) 


with 7^ the constant-curvature (/C = 0, ±1) spatial 
metric, let W* = dx M /dA, with A an affine param¬ 
eter, be a lightlike tangent to a geodesic curve, i.e. 
U fl U„ = 7 ^a 2 U l W - (U 0 ) 2 = 0 and = 

0. Expanding the geodesic equation in terms of the 
connections associated with the metric (1) and tak¬ 
ing into account the lightlike character of U , one 
finds that 


d LA 0 
~d\ 


+ H(U 0 ) 2 = 0, 
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which implies that 


stein equations read 


d 

dA 




= 0 , 


( 2 ) 


where the Hubble scale is H = d In a/di. Eq. (2) is 
solved by choosing the affine parameter A to satisfy 
dA = [a(t)/a(t{)]dt, with if a reference time, today 
say. We now assume our spacetime to begin at some 
initial coordinate time fj, which can take any value 
between 0 say, to —oo; this depends on the actual 
cosmological realization. The average Hubble rate 
along the geodesic parameterized by A is found to 
be 

s m - m A ( „, - ff(A)dA 

1 f 1 _ a[Hti)} \ 

A(if) — A(fi) \ a [A(if)] J 

<- (3) 

-A(i f )-A (i;)’ W 


so that in order for H ave ra ge to be strictly positive, 
a condition which is generally satisfied in inflation¬ 
ary models, one finds that the interval in affine pa¬ 
rameter must be finite, and therefore that the space- 
time under consideration is not geodesically com¬ 
plete. This argument cane be extended to timelike 
geodesics and more arbitrary cosmological models, 
i.e. with no specific assumptions regarding homo¬ 
geneity and isotropy. This requires the definition of 
a local expansion rate that is not dependent on the 
special FL metric solution; in this case, it is the de¬ 
viation between neighboring geodesics that needs be 
used explicitly to define the expansion rate (in the 
highly symmetric FL universe, the geodesic devia¬ 
tion is given by the expansion only, as this is the only 
relevant observable). The conclusion then is that, 
regardless of any energy condition, inflating space- 
times are past incomplete. 

An obvious way out of this problem consists in al¬ 
lowing the average Hubble rate to be negative. This 
amounts to having some amount of contraction, and 
hence, given that we observe the Universe to be cur¬ 
rently expanding, that it must have bounced. In the 
framework of GR however, this is not always easy. 

Using the metric (1) and a fluid stress-energy ten¬ 
sor = (p + P) u^u v + Pg^ with energy den¬ 
sity p , pressure P, and u ; , a timelike vector, the Ein- 


H 2 H —2 — oP’ H + H 2 — - — — -(p + 3P), 
a z 3 a o 

(4) 

leading to 

h = 54 ( ' ,+p >- <5) 

(we use natural units where h = c = 8i rG N = 1 
so that the Planck mass Mpi = G -1 / 2 is dimension- 
less). Although having an inflationary phase with 
a > 0 merely demands the violation of the Strong 
Energy Condition (SEC: p + 3P > 0), a bounce, 
requiring H —>■ 0 while H > 0 at the same time, 
implies that either the spatial sections must be pos¬ 
itively curved (/C > 0) or the Null Energy Condi¬ 
tion (NEC: p + P > 0) must be violated. In the 
former case, the scale factor at the bounce a B is ob¬ 
tained as the solution of 3/C/a 2 = p(a B ) and must 
satisfy /C/a 2 > — P(a B ). This condition is for in¬ 
stance fullfilled in the very simple case in which a 
single scalar field evolves in a potential with a local 
maximum. [10, 11], The bouncing solution, seen 
generically as indicated on Fig. 1, does have an ac¬ 
celerating phase, the scale-factor curve being convex 
at the bounce; although this technically implies the 
SEC to be also violated during a bouncing epoch, 
this cannot be understood as an inflating phase since 
the accelerating phase is not associated with a large 
increase of the size of the Universe. 

In the more familiar (to inflation-oriented cosmol- 
ogists) case of vanishing or negligible spatial curva¬ 
ture 1 , as mentioned above, Eqs. (4) imply a much 
more stringent constraint, namely that the NEC be 
violated; as discussed in Sec. 4., this often leads to 
various instabilities which then need to be tamed in 
order for the model to make any sense at all. 

3. Standard model puzzles, bouncing solutions 
- new issues 

The reason why the inflationary scenario is so fash¬ 
ionable stems from its successes in solving the stan¬ 
dard hot big bang puzzles in a unified way, while 

’Note that this is an assumption that can only be checked a 
posteriori: given a material content with positive and negative en¬ 
ergy components, one must first solve the Friedmann equation for 
the minimum scale factor a B , and then verify that the curvature 
term /C/a| is indeed negligible with respect to all other contri¬ 
butions. 
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Time t 


Figure 1: Typical time evolution of the scale factor a(t) 
(dashed line), Hubble rate H(t) = a/a (full line) and Hubble 
length £ h ( t ) = H ~ 1 (dotted line) for a bouncing scenario. For 
large negative times (conventionally setting the bounce at t = 0), 
the scale factor decreases in a non-accelerated way, then it curves 
up, accelerating and rendering the curve convex, finally con¬ 
necting, not necessarily in a symmetric way, to a more standard 
non-accelerating expansion. The Hubble rate starts vanishingly 
small, then decreases to large negative values, passes through the 
bounce almost linearly increasing, reaches a maximum and then 
decreases back to its usual behavior. The Hubble length is orig¬ 
inally very large, reaches a minimum and diverges at the bounce 
point: there is no super-Hubble scale at the bounce! 


at the same time providing a means of producing 
perturbations whose spectrum can be made to agree 
with all known data. Can a bouncing scenario, on 
top of naturally avoiding the singularity, propose sat¬ 
isfying solutions to the standard hot big bang puz¬ 
zles? If it is the case, can a bouncing scenario pro¬ 
vide a means to generate cosmological perturbations 
whose statistics agree with observations? As men¬ 
tioned before, we refer the reader to ref. [7] for a de¬ 
tailed discussion of this latter question, and focus in 
the remainder of this review on bouncing solutions 
to the background cosmological problems and on a 
review of existing bouncing models. 

3.1. Horizon and flatness puzzles 

The standard Hot Big-bang model suffers from 
a few puzzling problems, and we will treat in this 
section how a bounce, which implies a contracting 
phase preceding the current expansion, deals with 


the two most important, namely the horizon and flat¬ 
ness problems. We refer the reader to Refs. [7, 12] 
for more details relative to the other commonly ad¬ 
dressed puzzles. 


— Horizon 


The horizon problem relates to the inability to ex¬ 
plain the quasi-homogeneity of the observable uni¬ 
verse within the context of standard cosmology in 
which the entire evolution of the universe consists in 
decelerated expansion during the radiation and mat¬ 
ter dominated epochs. In the context of inflation, the 
necessity of a period of accelerated expansion for a 
period lasting a minimum of N ~ 60 e-folds can be 
phenomenologically understood by computing the 
solid angle subtended by causally connected regions. 
Assuming an expansion history in which the uni¬ 
verse is initially radiation-dominated, then dominated 
by the inflaton, represented by a fluid X with equa¬ 
tion of state parameter w x until a redshift z e , then 
once again radiation-dominated until the last scat¬ 
tering surface at Z| ss and finally matter dominated 
till today, we find that [ 13] 


-1/2I 1 


1 — (1 + ~lss) ^ 

1 — 3 w x 1 


L (l + ^ls S ) V2 11 


where N = ln(o e /oi) with a; the scale factor at the 
onset of the A - -dominated period. If we assume that 
N = 0, we recover standard cosmology and Af l ~ 
0.85 degrees. This corresponds to a total of about 
10 6 causally disconnected regions in which, weirdly 
enough, the CMB is everywhere the same up to 1 
part in 10 5 . Increasing Afl is easily achieved in the 
context of inflation by requiring w x < —1/3 and 
large positive N. 

A similar calculation can be done in the context 
of bouncing cosmology. Here again, we shall as¬ 
sume a phase dominated by a fluid with equation of 
state parameter w x during which the Universe first 
contracts. The bounce is assumed non-singular, oc¬ 
curring at a redshift z\ } and short enough that we can 
ignore its contribution. It is followed by the stan¬ 
dard radiation and matter dominated phases. The 
solid angle subtended by causally connected regions 


4 
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is then 


Aft = 


1 

2 L 


1 — (1 + ^lss) ^ (1 + Zlss) ^ 


1 + 


1 + +ss 3(1 + W x ) c -N(l+3w x )/2 
1 + Zb [ 1 + 3w x 


2 (2 + 3w x ) 
1 + 3 w x 



(7) 


In this expression, N = ln(ai/a B ) < 0 with cq and 
a B the scale factors at the onset of the A-dominated 
contraction and at the bounce respectively. Here, 
large values of AO are obtained for w > —1/3 and 
large values of \N\. Thus, in contrast with inflation, 
the horizon problem can be solved using a fluid that 
satisfies all energy conditions. 


— Flatness 


The flatness problem is easily understood by work¬ 
ing with the ratio of the Friedmann equation with 
the critical density p cr ; t = 3 H 2 , which is the energy 
density the Universe would have if it had exactly flat 
spatial sections. In the presence of the spatial curva¬ 
ture term, the Friedmann equation (4) takes the form 

a 2jj2 = ^2 “ 1 = - 1 > ( 3 ) 


with 



Pcrit 


H 2 (t 0 ) 

H 2 (t) 


O.j(fo) 


-3(l+!«i) 


(9) 


From observations, we know that the total density 
parameter today is fl T (fo) — 1. It is easy to recast 
Eq. (8) in the convenient form [13] 

N / s -3(l+!«j) 


(a) =$>,(£ o) - 


3 =1 
N 


a 0 




V«o 




a 

V«o 


-3(1+1/+ 


-2 n -1 


( 10 ) 


At early times (small scale factor), the Universe is 
radiation dominated, and Eq. (10) simplifies to 

2 


n T (t) - i ~ 


Or (to) - 1 / 1 


Hrad(fo) \1 + Z 


(ID 


For z 1, f2(f) — 1 must be much less than 1. 
For instance, taking z nuc i = 3 x 10 s , U rad (£ 0 ) = 
10~ 4 , and fl T (f 0 ) — 1 = 0.01, one finds fl T (f nU ci) — 
1 ~ 10~ 15 . The value of the total density parame¬ 
ter at nucleosynthesis required to satisfy today’s ob¬ 
served value fl T (fo) ~ 1 is highly fine-tuned and 
thus highly improbable. This embodies the flatness 
problem of standard cosmology. 


Let us now consider the case of a bouncing uni¬ 
verse which contracts in a phase dominated by a fluid 
of equation of state parameter w x , bounces and then 
expands according to the standard scenario. Note 
that the set of equations above do not apply at the 
bounce point where H = 0. In fact, in the pres¬ 
ence of a spatial curvature term, fl T diverges at the 
bounce. 

For a universe dominated by some fluid X, one 
has 


tV fl x (fi) — [!7 T (fi) — 1] (a/a{) 1+3wx ' 

( 12 ) 

The total density parameter at the end of the con¬ 
tracting phase at t~ is given by 


fi T (i_) - 1 ~ 


n T (fr)-i 

fix(fi) 


_ \ l+3«ix 


(13) 


while it is given by Eq. (11) at the beginning of the 
expanding phase, for t = t + and z = z + . The dif¬ 
ference 


Afi T = [n T (* + )-i]-[fi T (t_)-i] (14) 


can be computed using Eq. (8) and the Taylor expan¬ 
sion of the scale factor close to the bounce. 



One finds 


AU t 


3 P (U\ 2 

2 UJ 


(15) 

(16) 


Generically, in the absence of any fine-tuning, one 
should assume f2 T (fi) — 1 and fi.Y(oi) take values 
of 0(1) while it is known that fl rad (fo) — 10 -4 and 
fl T (fo) — 1 < 10~ 2 . Hence we have 


1+3 wx 


3/3 

2 



< 10 6 x . (17) 
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Taking as before z + ~ 10 28 , and with w = 1/3, we 
have 

^ < HT 50 , (18) 

where N < 0. Thus, for N < —60, and j3 of order 
1, a bounce with a short characteristic timescale and 
a large value of the scale factor at the bounce such 
that t c /a B < 10 -25 can satisfy current constraints 
on the spatial curvature of the universe. 

3.2. Shear/BKL instability 

In a bouncing scenario, the standard puzzles find 
natural solutions because of the contracting phase. 
However, such a phase can also induce another prob¬ 
lem: the fate of any initial amount of anisotropy. To 
focus on this question, we consider a spatially flat 
model whose dynamics derives from the Bianchi I 
metric ds 2 = —df 2 + a 2 {t)dx 2 , whose spatial part 
reads 

da: 2 = e 20 *^da: 2 + e^^d y 2 + e 2e *^d z 2 , (19) 


g 2AT 3/1 / t c 
2 \a B 


Friedmann equation reads 

rr2 _ 1 3/C PmO PrO Pd 0 P<j >0 

3 a 2 a 3 a 4 a 6 a 3(i+ui*) 

( 22 ) 

If this so-called ekpyrotic phase [17], lasts long en¬ 
ough, it eventually comes to dominate over all other 
constituents when a —> 0, including the shear con¬ 
tribution. The Universe then bounces and starts ex¬ 
panding again while in a fully symmetric FL phase, 
a condition absolutely required to explain the obser¬ 
vational data. 

The bounce itself is another matter, which we now 
turn to. 

4. Existing models 

There exists a large number of bouncing cosmo¬ 
logical models in the literature; we shall refer the 
reader to Ref. [7] for an exhaustive review and all 
the relevant references. We will instead focus here 
on a few models and give concrete examples. 


with = Qx + 0 y + 0 z = 0. Plugging (19) into 
the Einstein equations generalizes (4) to 

+ ( 20 ) 

x ' i 

and 

H = -\{p + P)-\^l (2D 

i 

where we have identified the shear energy density pg 
contained in the anisotropy stemming from the func¬ 
tions dp. Eqs. (20) and (21) imply that 0i + 3H9i = 
0, and therefore pg ex a -6 . 

With dust and radiation scaling as p m oc a~ 3 and 
p r ex a~ 4 respectively, the above result is a catastro¬ 
phe: as the universe contracts, any initial anisotropy, 
however small 2 , will grow until it eventually dom¬ 
inates the dynamics. This was shown [16] by Be¬ 
linsky, Khalatnikov and Lifshitz (BKL) to induce an 
instability sufficient to spoil the bounce. 

One way out of the shear problem is to add an 
extra component, usually a scalar field in a poten¬ 
tial satisfying specific constraints, with large effec¬ 
tive equation of state 1, so that the resulting 

2 Actually, the problem only arises in the presence of primor¬ 
dial classical shear: it has been shown that if the primordial shear 
is generated by quantum vacuum fluctuations, scalar and vector 
perturbations remain comparable [14, 15]. 


4.1. Classical bounces 

Bouncing models predate by many decades the in¬ 
flationary paradigm, as they were first introduced in 
the 1930’s. Classical models involve unconventional 
perfect fluids or scalar fields with possibly non stan¬ 
dard kinetic terms, or various combinations of those. 
The most conservative setup that may be used to ob¬ 
tain a bounce is to introduce spatial curvature and to 
violate the strong energy condition. In such a setup, 
scalar field matter is required in order to achieve 
p + 3 P < 0 and either a quasi-symmetric bounce 
or a phase of inflation is needed to drive towards 
zero. 

• Perfect fluids 

In a theory restricted to GR and FL spacetime, gener¬ 
ic ally, for /C +1, the null energy condition has to 
be violated in order to obtain a bounce, as discussed 
below Eq. (5). Exotic hydrodynamical fluids that vi¬ 
olate the null energy condition, and thus all other 
energy conditions, are a priori allowed, and mod¬ 
els using those can be built in the framework of an 
FL spacetime. Assuming an expansion for the scale 
factor of the form 

a = a 0 + b V 2n + dif n+1 + er ] 2n+2 , (23) 
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where n > 1, and (cio, b, d , e) constant, it is possi¬ 
ble to compute, in a fully analytic way, the evolution 
of adiabatic perturbations around the bounce [18], 
The possible choices for n and 1C are: (i) n > 1 and 
1C 0; (ii) n > 1 and K. = 0; (iii) n = 1 and 
d ^ 0 V/C; (iv) n = 1 and d = 0 V/C; setting d = 0 
restricts to a symmetric bounce. In the first three 
cases, the Bardeen potential $ describing the gauge- 
invariant perturbations turns out to be singular at the 
bounce, while in case (iv), although >1> is well be¬ 
haved, the NEC needs be violated even if /C = +1. 
Given that at late times, it has to be satisfied, there 
must exist a time i* at which p(t*) + P(f*) = 0. 
It turns out that, at this NEC transition, the growth 
of *I> is unlimited, raising potential questions on the 
perturbative expansion through a bouncing phase. 

When entropy perturbations are considered in ad¬ 
dition to the adiabatic ones, they are found to be 
sourced by the interaction of the hydrodynamical flu¬ 
ids involved in the cosmological evolution. This im¬ 
plies that the fluids do not evolve independently. In¬ 
clusion of entropy perturbations has the effect of reg¬ 
ularizing the Bardeen potential and its derivatives at 
the NEC transition. It may thus be concluded that 
perfect fluid dominated bouncing models in which 
both adiabatic and entropy perturbations are taken 
into account are regular and do not necessarily lead 
to strong backreaction effects of the perturbations 
onto the background geometry [19], 

When it comes to violating energy conditions, it 
is tempting to make use of scalar fields: most imple¬ 
mentations of the inflationary paradigm are based on 
scalar fields, and the bounce does no better in that re¬ 
spect! Bouncing models powered by such fields can 
be broadly distinguished in two categories, depend¬ 
ing on the coupling (minimal or extended) with the 
geometry. We now discuss both these possibilities. 

• Minimally coupled scalar fields theories 

Theories using minimally coupled scalar fields can 
be separated in two categories. First, an ordinary 
scalar field, with standard kinetic term and a po¬ 
tential. In this case, in order to obtain a bouncing 
cosmology and preserve the weak energy condition, 
one needs 1C = + 1. We shall not dwell here with 
such cases, which then demand the curvature prob¬ 


lem to be addressed independently, and may result in 
perturbations being large and potentially highly non 
Gaussian [20], 

The second category naturally involves non stan¬ 
dard kinetic terms; those can be generalized to the 
Galileon theories (non minimal coupling, see below). 
Their advantage over the standard kinetic terms is 
that those can be implemented in a flat FL universe. 

— Ghost condensates 

The simplest possible example of a non standard ki¬ 
netic terms consists in merely switching its sign, mak¬ 
ing it a so-called ghost field, which yields an un¬ 
tenable theory because instabilities, both classical 
and quantum, will immediately develop and ruin any 
configuration. Although one can use a ghost as an ef¬ 
fective means to initiate a bounce with 1C = 0 [21], 
it makes more sense to induce NEC violations with 
dynamical ghosts. This can be achieved in higher- 
derivative theories with second order equations of 
motion that by construction prevent gradient or ghost 
instabilities. This is the ghost condensate mecha¬ 
nism [22], whose features we can sketch with the 
Lagrangian 

C = P{X) where A' = (24) 

and where the pressure P is an arbitrary function of 
the kinetic energy X. Eq. (24) in a flat FL metric 
yields 

j t ( a3 P,x& =0, (25) 

where P x = d/'/dA. 

If A is a constant and P x = 0 at X = X c , the 
equation of motion yields the solution 

<t>=y/2X c t. (26) 

Given that 

p + P = 2XP x , (27) 

and the constraint X > 0, a violation of the NEC 
can take place if P x < 0 in some interval of the 
values of A'. This, the condition that P x = 0 at 
X = X c and requiring that P xx > 0 at X = X c in 
order to prevent the existence of ghosts implies that 
the function P(X) should have a local minimum at 
X c , (see Fig. 2). This construction was employed 
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in the so-called K —bounce [23], new Ekpyrotic [24] 
and matter bounce [25] scenarios. Unfortunately, the 
ghost condensate phase in models such as this can¬ 
not be smoothly connected with a branch P x > 0 
at X = 0. Ghost condensate type models there¬ 
fore do not admit a stable Poincare-invariant vacuum 
state and are thus severely flawed. Consequences of 
this instability were for instance demonstrated ex¬ 
plicitly in the case of the new Ekpyrotic scenario in 
Refs. [26, 27], 

— Ekpyrotic potential and ghost condensation 

As was mentioned earlier, the contracting phase of 
a bouncing cosmology is generically unstable un¬ 
der the growth of anisotropies and leads to chaotic 
mixmaster oscillations unless a period of ekpyrotic 
contraction with w > 1 is invoked or if the con¬ 
traction is sufficiently brief. The smooth transition 
from ekpyrotic contraction to expansion through a 
non-singular bouncing phase relying on the NEC vi¬ 
olation (with w < 1) was first studied in the new 
Ekpyrotic model [24], It is obtained with a single 
scalar field rolling down a steep negative potential 
during the ekpyrotic phase and then undergoing a 
ghost condensation. In this approach the ghost con¬ 
densate Lagrangian is thus supplemented with a po¬ 
tential term V(cf>). The function P(X ) realizing the 
ghost condensate phase, and the ekpyrotic potential 
V(o) are depicted in Fig. 2 As shown for instance 
in [28] however, this model suffers from a gradi¬ 
ent instability and from the regrowth of the initial 
anisotropy during the bouncing phase. In addition, 
the absence of a Lorentz-invariant vacuum remains, 
as in the ghost-condensate model. It also predicts a 
blue spectrum of curvature perturbations. 

• Conformal Galileons 

Instead of realizing NEC-violations by relying on 
Lagrangians that are restricted to general functions 
of the kinetic term only, it is possible to construct 
yet more general NEC-violating theories with La¬ 
grangians that exhibit couplings between various sca¬ 
lar field derivative terms [29]. Such scalar field the¬ 
ories are called Galileon theories and have drawn 
much interest due to the fact that they naturally ad¬ 
mit self-accelerating solutions. In Galileon theories, 



Figure 2: Top: ghost condensate kinetic function P(X). Bot¬ 
tom: ekpyrotic potential Plots obtained from Ref. [28], 


the scalar Lagrangian involves higher derivative in¬ 
teractions with at most second order derivatives in 
the equations of motion and is invariant under (pos¬ 
sibly conformal or DBI-conformal) Galilean trans¬ 
formations. Galileon theories are a subclass of the 
theory of Generalized Galileons which are described 
by the most general scalar-tensor (Horndeski) action 
leading to second order equations of motion. Denot¬ 
ing the scalar field by 7r(a"), it reads 


Sn=J d 4 Xy/^g^P{n,X) + Ga(x,X)an 
+ G K {ir, X)K + Gn, x [( n ^) 2 - (V M V„7 t ) 5 


+ G g (ir,X) G^WW'i r- 


G g 

6 


(□t r) 3 - 3D7T (V M V„7r) 2 + 2 (V^V^tt) 5 


(28) 


with 1Z the scalar curvature, the functions Go, G-ji 
and Gg being arbitrary functions of the field n and 
its kinetic energy X ; the last two functions make 
manifest the non minimal couplings with the grav- 
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itational sector. 

Conformal formulations of Galileon theories are 
particularly advantageous because the 4D conformal 
group reduces to the Anti de Sitter Minkowski and 
de Sitter symmetry groups for particular solutions of 
the dilaton equation of motion [30]. More impor¬ 
tantly, whereas superluminal propagation of pertur¬ 
bations is common in higher derivative theories such 
as Galileons, in conformal Galileon theories [30], 
perturbations of the scalar can be shown to travel 
with a speed at most equal to the speed of light in 
the entire phase space as long as matter fields are 
excluded [31]. In [30] a simple and almost viable 
example is provided, with the theory described by 
the Lagrangian 


C = / 2 e 27r (S7r) 2 +^3(57r) 2 D7r+^3-(l-|-a)(57r) 4 , 

(29) 

where /, A and a are constant. Provided C features 
a negative kinetic term, as is the case in Eq. (29), this 
theory admits a time-dependent de Sitter solution. 


e 


7T 


1 

-Hot' 


2 1 A 3 

with H 0 = ~— -- (30) 

3 1 + a j 


with A the theory’s strong coupling scale and vio¬ 
lates the NEC. The sound speed is subluminal for 
0 < a < 3 but the NEC violating solution, as is 
the case for the ghost condensate cannot be smoothly 
connected with a Lorentz-invariant vacuum solution. 
This theoretical setup nevertheless does bare a par¬ 
ticularly interesting aspect which was dubbed the 
“Galileon Genesis”, namely that the cosmological 
solution displayed here is an attractor solution. This 
allows the possibility to have either emerging cos¬ 
mological evolution or bouncing solutions. 

4.2. Semi-classical and quantum bounces 

Semi-classical models are those involving quan¬ 
tized scalar fields in classical spacetimes. The vac¬ 
uum state being ill defined on a curved background, 
except in the adiabatic limit of slowly varying scale 
factor, a regularization or renormalization scheme is 
required to cure such semi-classical theories of the 
infinities that arise in the formal expression for the 
stress-energy tensor. These infinities are associated 
with the inability to properly define the creation and 
annihilation operators and thereby unambiguously 
remove the infinite vacuum term in the usual way. 


Renormalization results in the inclusion of higher or¬ 
der curvature (counter-) terms in the Einstein-Hilbert 
action leading to new terms in the Friedmann equa¬ 
tions and to the possibility of constructing singular¬ 
ity avoiding cosmologies. 

A more ambitious approach to describe quantum 
gravitational effects is string theory. The full ac¬ 
tion of superstring theory possesses scale factor du¬ 
ality and time reversal symmetry, which can be used 
to construct non-singular cosmologies. These re¬ 
quire a branch change that smoothly interpolates be¬ 
tween contracting and expanding spacetimes. This 
is the well-known pre-big-bang cosmology ’ [32], In 
its original version this model consists of only the 
dilaton field and the metric. At tree level, it can 
be shown that a contracting cosmology in the string 
frame corresponds to an expanding cosmology in the 
Einstein frame. These two frames are simply related 
by a conformal transformation. It is therefore not 
clear whether to identify such a tree-level cosmol¬ 
ogy with a bouncing cosmology per se. However, 
with the inclusion of loop corrections, it is possible 
to show that the cosmology is indeed non-singular, 
and it is then plausible that one identifies this non¬ 
singular evolution with branch changing. 

Another possibility to smooth out the curvature 
singularity is the inclusion of a coupling to a mat¬ 
ter or radiation fluid in the tree level effective 4D 
action of string theory [33, 34]. In order to study the 
propagation of cosmological perturbations through 
a bounce, it can either be modeled by a discontinu¬ 
ity across a spacelike hypersurface, in which case it 
is singular, and the behavior of cosmological pertur¬ 
bations transfered through the bouncing phase de¬ 
pends on how the Israel junction conditions are im¬ 
plemented, or, alternatively, if one smooths the cur¬ 
vature singularity by including higher order correc¬ 
tions in the effective action, then it becomes possible 
to actually follow perturbations through the bounce. 

Another string way to a non singular cosmology 
involves the motion and interactions of higher di¬ 
mensional (mem)branes in yet higher dimensional 
bulk geometries. While initially, models were based 
on branes embedded in, e.g. 5-dimensional bulk ge¬ 
ometries, more recent models, based on either het¬ 
erotic M-theory or on the compactification of 10- 

3 At low energy, the four dimensional effective action of the 

ekpyrotic model is equivalent to a modified version of the pre-big 
bang model. 
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dimensional type II A/B superstring theory mani¬ 
folds, and on the stabilization of the “moduli” fields 4 , 
have led to interesting brane dynamics in warped 
parts of the geometry. In general, these construc¬ 
tions lead to additional terms in the Friedmann equa¬ 
tions or unconventional kinetic terms for the inflaton 
field that can lead to a period of inflation, to bounc¬ 
ing branes and to cyclic cosmologies. The ekpy- 
rotic model [35] is one realization of brane cosmol¬ 
ogy based on heterotic M-theory while examples of 
bouncing cosmologies in warped string compactifi- 
cations can be found in e.g. [36, 37]. 

A final option, also fully quantum, consists in as¬ 
suming the energy scale at which the bounce takes 
place to be sufficiently small that a Wheeler de Witt 
treatment of quantum cosmology would be appro¬ 
priate. Although one naturally faces measurement 
questions in such a context, there exists ways to treat 
both background and perturbations on an equal foot¬ 
ing; a full set of predictions to compare with current 
data is however not yet available since only toy mod¬ 
els have been written down, but it would seem that 
a consistent model should be attainable in the near 
future (see Ref. [38] and references therein). 

5. Conclusions 

Although inflation appears to largely dominate the 
field of primordial cosmology, due, in particular, to 
the fact that many implementations have predicted 
consequences rather similar to the presently avail¬ 
able observations, bouncing alternatives are not en¬ 
tirely ruled out. In addition, a contracting phase fol¬ 
lowed by a bounce can solve the primordial singu¬ 
larity problem, which renders such models attractive 
and worth investigating. It must be conceded how¬ 
ever that, currently, most of bouncing models have 
difficulties, either because they demand a compli¬ 
cated theoretical framework or because their predic¬ 
tions disagree with cosmological data. 

From a purely theoretical standpoint, one might 
argue that bouncing cosmologies relying on either 
ghost condensates or unknown quantum gravitational 
effects [7] in order to successfully avoid the clas¬ 
sical singularity should, in view of Occam’s razor, 
be disfavored when compared to much simpler in¬ 

4 Moduli fields are fields that appear after compactification. 
They are identifiable with the (a priori unfixed) sizes and shapes 
of cycles in the higher dimensional Calabi-Yau manifold which 
the theory lives on. 


flationary models involving perfectly well-behaved 
scalar fields [4]. One should nevertheless remember 
that the singularity problem of big-bang cosmology 
will need to be addressed at some stage, and its so¬ 
lution, however contrived it may look from our per¬ 
spective, may end up being quite natural within a 
few decades. In other words, one should not refute 
a theory or a paradigm on philosophical grounds but 
instead on whether it is able to answer as of yet unan¬ 
swered physical questions and whether it agrees with 
the data or not. 

From an observational perspective, given increas¬ 
ingly stringent constraints imposed by present day 
cosmological data, many bouncing models are un¬ 
der pressure as they naturally predict either exactly 
scale-invariant scalar perturbations, or even slightly 
blue spectra. Reddening the spectrum often demands 
new components which may contribute in a non neg¬ 
ligible way and produce unobserved isocurvature mo¬ 
des. Furthermore, the bounce itself, involving either 
NEC violating fields or positive spatial curvature, 
might induce large non Gaussianities [20], Present 
day bouncing models scarcely agree with all exist¬ 
ing observational constraints but it must be noted 
that neither do most inflationary models [4] . More 
work is needed to reach definite conclusions, per¬ 
haps along the lines of purely quantum models? 
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